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Abstract 

A unified treatment of mass varying dark matter coupled to cosmon-Zifce dark energy is shown 
to result in effective generalized Chaplygin gas (GCG) scenarios. The mass varying mechanism 
is treated as a cosmon field inherent effect. Coupling dark matter with dark energy allows for 
reproducing the conditions for the present cosmic acceleration and for recovering the stability 
resulted from a positive squared speed of sound c 2 s , as in the GCG scenario. The scalar field 
mediates the nontrivial coupling between the dark matter sector and the sector responsible for 
the accelerated expansion of the universe. The equation of state of perturbations is the same as 
that of the background cosmology so that all the effective results from the GCG paradigm are 
maintained. Our results suggest the mass varying mechanism, when obtained from an exactly 
soluble field theory, as the right responsible for the stability issue and for the cosmic acceleration 
of the universe. 
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I. INTRODUCTION 



The ultimate nature of the dark sector of the universe is the most relevant issue related 
with the negative pressure component required to understand why and how the universe 
is undergoing a period of accelerated expansion jl-5|. A natural and simplistic explana- 
tion for this is obtained in terms of a tiny positive cosmological constant introduced in the 
Einstein's equation for the universe. Since the cosmological constant has a magnitude com- 
pletely different from that predicted by theoretical arguments, and it is often confronted 
with conceptual problems, physicists have been compelled to consider other explanations 



for that 



121. 



Motivated by the high energy physics, an alternative for obtaining a negative pressure 
equation of state considers that the dark energy can be attributed to the dynamics of a scalar 
field 4> which realizes the present cosmic acceleration by evolving slowly down its potential 



V(4>) 



13 



14(. These models assume that the vacuum energy can vary [15j|. Following 



theoretical as well as phenomenological arguments, several possibilities have been proposed, 



such as fc-essence 



131, p 



rantom energy 



of modifications of gravity 
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191 ] , cosmon fields [20] , and also several types 



23]. 
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27] coupled to 



One of the most challenging proposals concerns mass varying particles 
the dark energy through a dynamical mass dependence on a light scalar field which drives 
the dark energy evolution in a kind of unified cosmological fluid. The idea in the well-known 
mass varying mechanism 12J, I26l-l28| is to introduce a coupling between a relic particle and 
the scalar field whose effective potential changes as a function of the relic particle density. 
This coupled fluid is either interpreted as dark energy plus neutrinos, or as dark energy 
plus dark matter 29|, |30(. Such theories can possess an adiabatic regime in which the scalar 



field always sits at the minimum of its effective potential, which is set by the local mass 
varying particle density. The relic particle mass is consequently generated from the vacuum 
expectation value of the scalar field and becomes linked to its dynamics by m(4>). 

A scenario which congregate dark energy and some kind of mass varying dark matter in a 
unified negative pressure fluid can explain the origin of the cosmic acceleration. In particu- 
lar, any background cosmological fluid with an effective behaviour as that of the generalized 
Chaplygin gas (GCG) 7H9J naturally offers this possibility. The GCG is particularly rel- 
evant in respect with other cosmological models as it is shown to be consistent with the 
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observational constraints from CMB 



35[, and gamma ray bursts 



311 ]. supernova 32l434j| . gravitational lensing surveys 



Moreover, it has been shown that the GCG model can be 
accommodated within the standard structure formation mechanism {7,9, 



33]. 



In the scope of finding a natural explanation for the cosmic acceleration and the corre- 
sponding adequation to stability conditions for a background cosmological fluid, our purpose 
is to demonstrate that the GCG just corresponds to an effective description of a coupled 
fluid composed by dark energy with equation of state given by p(<fr) = —p{4>) and a cold dark 
matter (CDM) with a dynamical mass driven by the scalar field <fi. Once one has consistently 
obtained the mass dependence on 0, which is model dependent, it can be noticed that the 
cosmological evolution of the composed fluid is governed by the same dynamics prescribed 
by cosmon field equations. It suggests that the mass varying mechanism embedded into the 
cosmon-Z^e dynamics reproduces the effective behaviour of the GCG. In addition, coupling 
dark matter with dark energy by means of a dynamical mass driven by such a scalar field 
allows for reproducing the conditions for the present cosmic acceleration and for recovering 
the stability prescribed by a positive squared speed of sound c 2 s . At least implicitly, it leads 
to the conclusion that the dynamical mass behaviour is the main agent of the stability issue 
and of the cosmic acceleration of the universe. 

At our approach, the dark matter is approximated by a degenerate fermion gas (DFG). 
In order to introduce the mass varying behaviour, we analyze the consequences of coupling 
it with and underlying dark energy scalar field driven by a cosmon- type equation of motion. 
We discuss all the relevant constraints on this in section II. In section III, we obtain the 
energy density and the equation of state for the unified fluid and compare them with the 
corresponding quantities for the GCG. In section IV, we discuss the stability issue and 
the accelerated expansion of the universe in the framework here proposed. The pertinent 
comparisons with a GCG scenario are evaluated. We draw our conclusions in section V by 
summarizing our findings and discussing their implications. 



II. MASS VARYING MECHANISM FOR A DFG COUPLED TO COSMON- LIKE 
SCALAR FIELDS 

To understand how the mass varying mechanism takes place for different particle species, 
it is convenient to describe the corresponding particle density, energy density and pressure 



3 



as functional of a statistical distribution. This counts the number of particles in a given 
region around a point of the phase space defined by the conjugate coordinates: momentum, 
/3, and position, x. The statistical distribution can be defined by a function /(g) in terms of 
a comoving variable, q — a where a is the scale factor (cosmological radius) for the flat 
FRW universe, for which the metrics is given by ds 2 = dt 2 — a 2 (t)S ij dx t dx j . In the flat FRW 
scenario, the corresponding particle density, energy density and pressure are thus given by 

1 



n(a) 

Pm{a-,4>) 
Pm{a>,4>) 



7r 2 a 4 

1 

3vr 2 a 4 



dqq 2 f{q), 
dqq 2 {q 2 + m 2 (0)a 2 ) l/2 f(q), 
dqq' (q 2 + m 2 (<f>)a 2 )- l/2 f( q ). 



(1) 



37|]. For 



where the last two can be depicted from the Einstein's energy-momentum tensor 
the case where /(g) is a Fermi-Dirac distribution function, it can be written as 

f( q ) = {exp[(q-q F )/T } + l}- 1 , 

where T is the relic particle background temperature at present. In the limit where T tends 
to 0, it becomes a step function that yields an elementary integral for the above equations, 
with the upper limit equal to the Fermi momentum here written as q F = a /3(a). It results in 
the equations for a DFG 38] . The equation of state can be expressed in terms of elementary 
functions of /3 = /3(a) and m = m(</>(a)), 



n(a) 



3tt 2 ^ 
1 



Pm(a) = — 
87T 2 



/3(2/3 2 + m 2 ) a//3 2 + m? - arc sinh (/3/m) 



(2) 



Pm{a) 



1 

hi 2 



j3{-f3 2 - m 2 ) ^Jf3 2 + m 2 + arc sinh (/3/m) 
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One can notice that the DFG approach is useful for parameterizing the transition between 
ultra-relativistic (UR) and non-relativistic (NR) thermodynamic regimes. It is not manda- 
tory for connecting the mass varying scenario with the GCG scenario. 

Simple mathematical manipulations allow one to easily demonstrate that 

dp m (a 



n(a) 



and 



m(a) 



dn(a) 

dpm(a) 
dm(a) 



{.Praia) +Pm{a)), 



(Pm(a) - 3p m (a)), 



(3) 



(4) 



Noticing that the explicit dependence of p m on a is intermediated by (3(a) and m(a) = 
m((j)(a)), one can take the derivative of the energy density with respect to time in order to 
obtain 

• dp m (a) . dp m (a) 
p m = /3(a) +m{a)- 



8(3 (a) dm{a) 

■ , dp m (a) . dp m (a) 
71(a) — h m(a)- 



dn(a) dm(a) 

a dp m {a) -dmdpmia) 
= -3-n(o) +0 — , 5 

a an(a) d0 om(a) 

where the overdot denotes differentiation with respect to time (' = d/dt). The substitution 

of Eqs. (EH]) into the above equation results in the energy conservation equation given by 

. dim 

p m + 3H(p m + Pm) ~ ^"^0 ~ 3 ^™) = ' ( 6 ) 

where H = a/ a is the expansion rate of the universe. If one performs the derivative with 
respect to a directly from p m in the form given in Eq. (j2J), the same result can be obtained. 
The coupling between relic particles and the scalar field as described by Eq. (jBJ) are effective 
just for NR fluids. Since the strength of the coupling is suppressed by the relativistic increase 
of pressure (p ~ 3p), as long as particles become relativistic (T(a) = T /a » m(4>(a))) 
the matter fluid and the scalar field fluid tend to decouple and evolve adiabatically. The 
mass varying mechanism expressed by Eq. translates the dependence of m on into a 
dynamical behaviour. In particular, for a DFG, the consistent analytical transition between 
UR and NR regimes and their effects on coupling dark matter (m) and dark energy (0) are 
evident from Eq. ([6]). The mass thus depends on the value of a slowly varying classical scalar 



field m 



40| which evolves like a cosmon field. The cosmon- type equation of motion for the 



scalar field <fi is given by 

+ 3F0+^M = Q(# (7) 

where, in the mass varying scenario, one identifies Q{4>) as —(dm/d<j))/(dpm/dm). The 
corresponding equation for energy conservation can be written as 

P,+m P ,+P,)+^^ = o. (8) 

which, when added to Eq. (jBj), results in the equation for a unified fluid (p,p) with a dark 
energy component and a mass varying dark matter component, 

p + 3H(p + p) = 0, (9) 



where p = p (f> + p m and p = p^ + p m . 

As we shall notice in the following, this unified fluid corresponds to an effective description 
of the universe parameterized by a GCG equation of state. 

III. DECOUPLING MASS VARYING DARK MATTER FROM THE EFFECTIVE 
GCG 



Irrespective of its origin, several studies yield convincing evidences that the GCG sce- 
nario is phenomenologically consistent with the accelerated expansion of the universe. This 



31| given by 



scenario is introduced by means of an exotic equation of state [7j, [9 

P =-A s (^y, (io) 

which can be obtained from a generalized Born-Infeld action Q. The constants A s and a 
are positive and < a < 1. Of course, a = corresponds to the ACDM model and we 
are assuming that the GCG model has an underlying scalar field, actually real 0, Q] or 
complex 8, |9[. The case a = 1 corresponds to the equation of state of the Chaplygin gas 
scenario 7] and is already ruled out by data [3l|- Notice that for A s = 0, GCG behaves 
always as matter whereas for A s = 1, it behaves always as a cosmological constant. Hence 
to use it as a unified candidate for dark matter and dark energy one has to exclude these 
two possibilities so that A s must lie in the range < A s < 1. 

Inserting the above equation of state into the unperturbed energy conservation Eq. 
one obtains through a straightforward integration 



rnpert 



P = Po 



A s + 



:i-a s ] 

a 3(i+a) 



and 



p = -A„p 



A s + 



A -A, 



l/(l + Q ) 



- -a/(l + a) 



(11) 



(12) 



a s(i+a) 

One of the most striking features of the GCG fluid is that its energy density interpolates 
between a dust dominated phase, p oc a~ 3 , in the past, and a de-Sitter phase, p = —p, at 
late times. This property makes the GCG model an interesting candidate for the unification 
of dark matter and dark energy. Indeed, it can be shown that the GCG model admits 
inhomogeneities and that, in particular, in the context of the Zeldovich approximation, 
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these evolve in a qualitatively similar way as in the ACDM model js|. Furthermore, this 
evolution is controlled by the model parameters, a and A s . 

Assuming the canonical parametrization of p and p in terms of a scalar field 4>, 

1 •. 



P 
P 



+ V, 

- v, 



(13) 



allows for obtaining the effective dependence of the scalar field <fi on the scale factor, a, and 
explicit expressions for p, p and V in terms of <f). 

Following Ref. [33], one can obtain through Eq. (jHJ) the field dependence on a, 

. -l -a/(a+l) 



2 (a) 



Pop- - A) 

/,3(a+l) 



A, 



a 



3(a+l) 



and assuming a flat evolving universe described by the Friedmann equation H 2 
H in units of H and p in units of p Crit = 3H 2 /87rG), one obtains 



4>(a) 



where it is assumed that 



3(a + l) 



In 



y/1 - A,(l - a^ +1 )) - VY^A 



(14) 
p (with 

(15) 



(ao = 1) 



1 



In 



l- vt^a s 



3(a + l) 

One then readily finds the scalar field potential, 

= -Aj+^p {[cosh (3(a + 1)0/2)] ^ + [cosh (3(a + 1)0/2)]"^ 



(16) 



(17) 

If one supposes that energy density, p, may be decomposed into a mass varying CDM compo- 
nent, p m , and a dark energy component, p^, connected by the scalar field equations (JT])-([HD, 
the equation of state (TlOl) is just assumed as an effective description of the cosmo logical 
background fluid of the universe. Since the CDM pressure, p m , is null, the dark energy 
component of pressure, p^, results in the GCG pressure, p = p§. Assuming that dark energy 
obeys a de-Sitter phase equation of state, that is, p${$) = — p</>(0), the dark energy density 
can be parameterized by a generic quintessence potential, p^) = U(<f>), since its kinetic 
component has to be null for a canonical formulation. It results in U(<f>) = —p^{4>) = P, 
where p is the GCG pressure given by Eq. ( fl2i) . By substituting the result of Eq.( Tl5|) into 
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the Eq.f ll2p . and observing that H 2 = p, with p given by Eq. ljTT]) . it is possible to rewrite 
the GCG pressure, p, in terms of <fi. It results in the following analytical expression for £/(</>), 



U{4>) = P(j) {4>) = -p^) 



A* cosh 



3(a + l) 



2a 
1 + a 



(18) 



2 

which is consistent with the result for V{4>) = (l/2)(p(4>) — p(4>)) from Eq. f JTTj) . Since 
Pcjti^) + £></>(</>) = 0, the Eq. §E§ is thus reduced to 

dU (0) dp m dm(0) _ Q 
d(/> 9m d0 

and the problem is then reduced to finding a relation between the scalar potential U(4>) 
and the variable mass m(cj>). From the above equation, the effective potential governing the 
evolution of the scalar field is naturally decomposed into a sum of two terms, one arising 
from the original quintessence potential U (</>), and other from the dynamical mass m(<f>). For 
appropriate choices of potentials and coupling functions satisfying Eq. ffT9l) . the competition 
between these terms leads to a minimum of the effective potential. For quasi-static regimes, 
it is possible to adiabatically track the position of this minimum, in a kind of stationary 
condition. The timescale for <ft to adjust itself to the dynamically modified minimum of 
the effective potential may be short compared to the timescale over which the background 
density is changing. In the adiabatic regime, the matter and scalar field are tightly coupled 
together and evolve as one effective fluid. At our approach, once we have assumed the 
dark energy equation of state as p$ = —pc/,, the stationary condition is a natural issue 
that emerges without any additional constraint on cosmon-type equations. In the GCG 
cosmological scenario, the effective fluid description is valid for the background cosmology 
and for linear perturbations. The equation of state of perturbations is the same as that of the 
background cosmology where all the effective results of the GCG paradigm are maintained. 
The Eq. ffTSl leads to p + p = p m + p m which, in the CDM limit, gives 

p(a) + p(a) = m(a) n(a) + p m {= 0) = — m(a) (5\a). (20) 

Since the dependence of m on a is exclusively intermediated by 0(a), i. e. m(a) = m(</>(a)), 
from Eqs. (TTT1) . ([12]) and ( Tl5l) . after some mathematical manipulations, one obtains 



m(4>) = m 



tanh (3(a + l)f) 
tanh (3(a + l)4a) 



2a 
1 + a 



(21) 
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which is consistent with Eq. (!T9|) once n oc a~ 3 . One can thus infer that the adequacy to 
the adiabatic regime is left to the mass varying mechanism which drives the cosmological 
evolution of the dark matter component. 

To give the correct impression of the time evolution of the abovementioned dynamical 
quantities driven by 0, in the Fig. [TJwe observe the behaviour of m(a) and U(a) in confront 
with 0(a) and V(a) of the GCG. In the Fig. [2] we verify how the energy density p and the 
corresponding equations of state u for the unified fluid, p m + P<t>, which imitates the GCG, 
deviates from the right GCG scenario. We assume that the mass varying dark matter behaves 
like a DFG in a relativistic regime (hot dark matter (HDM)) and in a non-relativistic regime 
(CDM). For mass varying CDM coupled with dark energy with = —pj>, the effective GCC 
leads to similar predictions for u, independently of the scale parameter a. The same is not 
true for HDM which, in the DFG approach, when weakly coupled with dark energy, leads 
to the same behavior of the GCG just for late time values of a (a ~ 1). 

As one can observe, the mass varying mechanism allows for reconstituting the GCG 
scenario in terms of non-exotic primitive entities: CDM and dark energy. Furthermore, as 
we shall notice, the dynamical mass expressed by Eq. ( 12TT) overpasses the problematic issue 
of stability. 

It is also important to emphasize that, as in the case of the Chaplygin gas, where a = 1, 
the GCG model admits a d-brane connection as its Lagrangian density corresponds to the 
Born-Infeld action plus some soft logarithmic corrections. Space-time is shown to evolve 
from a phase that is initially dominated, in the absence of other degrees of freedom on the 
brane, by non-relativistic matter to a phase that is asymptotically De Sitter. The Chaplygin 
gas reproduces such a behaviour. In this context, the explicit dependence of the mass on a 
scalar field is relevant in suggesting that the mass varying mechanism, when obtained from 
an exactly soluble field theory, can be the right responsible for the cosmological dynamics. 

IV. STABILITY AND ACCELERATED EXPANSION 



Adiabatic instabilities in cosmological scenarios was predicted [41( in a context of a mass 
varying neutrino (MaVaN) model of dark energy. The dynamical dark energy, in this ap- 
proach, is obtained by coupling a light scalar field to neutrinos but not to dark matter. Their 
consequent effects have been extensively discussed in the context of mass varying neutrinos, 
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in which the light mass of the neutrino and the recent accelerative era are twinned together 
through a scalar field coupling. In the adiabatic regime, these models faces catastrophic 
instabilities on small scales characterized by a negative squared speed of sound for the ef- 
fectively coupled fluid. Starting with a uniform fluid, such instabilities would give rise to 
exponential growth of small perturbations. The natural interpretation of this is that the Uni- 
verse becomes inhomogeneous with neutrino overdensities subject to nonlinear fluctuations 



42] which eventually collapses into compact localized regions. 

In opposition, in the usual treatment where dark matter are just coupled to dark energy, 
cosmic expansion together with the gravitational drag due to CDM have a major impact on 
the stability of the cosmological background fluid. Usually, for a general fluid for which we 
know the equation of state, the dominant effect on the sound speed squared c 2 arises from 
the dark sector component and not by the neutrino component. 

For the models where the stationary condition (cf. Eq. ( [19]) ) implies a cosmological 
constant type equation of state, p^ = — p^, one obtains c 2 s = — 1 from the very start of the 
analysis. The effective GCG is free from this inconsistency. The coupling of the dark energy 
component with dynamical dark matter is responsible for removing such inconsistency by 
setting c 2 s ~ ^ > 0. The exact behavior for dark energy plus mass varying dark matter fluid 
in correspondence with the GCG is exhibited in the Fig. [3] for different GCG a parameters. 
A previous analysis of the stability conditions for the GCG in terms of the squared speed 
of sound was introduced in Ref. 33J], from which positive c 2 s implies that < a < 1. These 
results are consistent with the accelerated expansion of the universe ruled by the dynamical 
mass of Eq. ( I2T1) which sets positive values for (1 + 3(p^ + Pm)/(p</, + p m )), as we can notice 
in the Fig. [3j 

For CDM (p << m) the unified fluid reproduces the GCG scenario. For HDM (p » m), 
in spite of not reproducing the GCG, the conditions for stability and cosmic acceleration 
are maintained. Fig. [3] shows that the GCG can indeed be interpreted as the effective result 
for the coupling between mass varying dark matter and a kind of scalar field dark energy 
which is cosmologically driven by a A-type equation of state, p^ = —p^. 
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V. CONCLUSIONS 



The dynamics of the cosmology of mass varying dark matter coupled with dark en- 
ergy dynamically driven by cosmon- type equations was studied without introducing spe- 
cific quintessence potentials, but assuming that the cosmological background unified fluid 
presents an effective behaviour similar to that of the GCG. 

We have comprehensively analyzed the stability characterized by a positive squared speed 
of sound and the cosmic acceleration conditions for such a dark matter coupled to dark energy 
fluid, that exists whenever such theories enter an adiabatic regime in which the scalar field 
faithfully tracks the minimum of the effective potential, and the coupling strength is strong 
compared to gravitational strength. The matter and scalar field are tightly coupled together 
and evolve as one effective fluid. The effective potential governing the evolution of the 
scalar field is decomposed into a sum of two terms, one arising from the original scalar field 
potential U (</>), and the other from the dynamical mass of the dark matter. 

The mass varying behaviour of the dark matter component was determined from the 
assumption of a kind of A-type dark energy dynamics embedded in an effective cosmological 
scenario which reproduces the cosmological effects of the GCG. It is equivalent to decoupling 
mass varying dark matter from the effective GCG concomitantly with assuming the dark 
energy equation of state as p^ = —p$. The adiabatic regime naturally occurs without any 
additional constraint on scalar field equations. The unified fluid description is valid for the 
background cosmology and for linear perturbations. The equation of state of perturbations 
is the same as that of the background cosmology where all the effective results of the GCG 
paradigm are maintained. 

Unfortunately, we cannot provide a sharp criterion on the potential and on the mass 
varying dependence on the scalar field to discriminate between these two possibilities: the 
GCG scenario or an effective unified fluid imitating the GCG via scalar fields driven by 
cosmon- type equations. Many results for specific quintessence potentials that are found in 
the literature, when they reproduce stability and cosmic acceleration, are recovered, and 
speculative predictions for new scenarios featuring other mass dependencies on scalar fields 
can be made. It remains open if the present approach can lead to a natural solution of 
the cosmological constant problem. In the meanwhile we take the cosmon model analogy 
as an interesting phenomenological approach, through which we can reproduce the main 
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characteristics of the GCG. 

Given the fundamental nature of the underlying physics behind the Chaplygin gas and 
its generalizations, it appears that it contains some of the key ingredients in the description 
of the Universe dynamics at early as well as late times. Our results suggest that the mass 
varying mechanism, when eventually derived from an exactly soluble field theory, which is 
noway trivial, can be the effective agent for the stability issue and for the cosmic acceleration 
of the universe, once it can effectively reproduce the main characteristics of a GCG scenario. 
It also stimulates our subsequent investigation of the evolution of density perturbations, 
instabilities and the structure formation in such scenarios. To summarize, we expect that 
the future precise data can provide more strong evidence to judge whether the dark energy 
is the cosmological constant and whether dark energy and dark matter can be unified into 
one cosmological background effective component. 
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FIG. 1: The cosmological evolution of m{a) and U(a) correlated with that of <p(a) and V{a) of the 
GCG as function of the scale factor (first plot) and the corresponding dependence on (f) (second 
plot). We are considering the GCG scenarios with A 3 = 3/4 and a = 1, 1/2, 1/4, 1/8. 
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FIG. 2: Energy densities, p, and equations of state, u = p/p, as function of the scale factor 
for the effective GCG fluid and its components: mass varying dark matter p m and cosmon-/i£;e 
dark energy pA. The results are compared with those of real GCG scenarios with A s = 3/4 and 
a = 1, 1/2, 1/4, 1/8. Dark matter is assumed to be described by a DFG in relativistic (p >> m) 
and non-relativistic regimes (p « m) at present. 
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FIG. 3: Squared speed of sound c 2 s = dp/dp and cosmic acceleration a/ (a p) as function of the scale 
factor for the unified fluid decomposed into mass varying dark matter and cosmon-Zi/ce dark energy 
components. Again the results are compared with those of real GCG scenarios for A s = 3/4 and 
a = 1, 1/2, 1/4, 1/8, for mass varying dark matter described by a DFG in relativistic (p >> m) 
and non-relativistic regimes (p « m) at present. 
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